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Direct photons have been measured in y?s NN — 200 GeV d+ Au collisions at midrapidity. A wide 
Pt range is covered by measurements of nearly-real virtual photons (1 < pr < 6 GeV/c) and real 
photons (5 < pr < 16 GeV/c). The invariant yield of the direct photons in d+Au collisions over 
the scaled p+p cross section is consistent with unity. Theoretical calculations assuming standard 
cold-nuclear-matter effects describe the data well for the entire pr range. This indicates that the 
large enhancement of direct photons observed in Au+Au collisions for 1.0 < pr < 2.5 GeV/c is due 
to a source other than the initial-state nuclear effects. 

PACS numbers: 25.75.Dw 



Direct photons in both Au+Au and p+p collisions were 
measured at the Relativistic Heavy Ion Collider [ll-Q over 
a wide pt range, which was achieved through measure- 
ments of both real photons and nearly-real virtual pho- 
tons [6]. For 1.0 < pt < 2.5 GcV/c, a significant ex- 
cess of direct photons over the binary-scaled p+p yield 
was observed in central Au+Au collisions, suggesting 
the existence of thermal photons emitted from the hot 
medium. The key to measurements of the direct pho- 
ton production for px < 5 GeV/c is the use of virtual 
photons, which greatly reduces the background of pho- 
tons from 7r°,?7 — > 27. For px > 4 GeV/c, real photons 
are used and previous Au+Au measurements [5j indi- 



cate agreement with the binary-scaled p+p collisions over 
4 < pt < 22 GeV/c. However, effects either in the initial 
state or in the medium created in Au+Au collisions may 
cancel, making the d+Au measurement crucial to under- 
standing the Au+Au results, because only initial-state 
effects are present in d+Au collisions. 

Cold-nuclear-matter (CNM) effects may play an im- 
portant role in direct photon production in A+A colli- 
sions and possibly modify the production rate compared 
to p+p collisions. CNM effects in the measured direct 
photon yield include interplays of various initial-state ef- 
fects such as the Cronin enhancement 0, isospin effect, 
modification of the nuclear parton distribution functions 
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(nPDFs) inside the nucleus 
ergy loss of colliding partons 




and the initial-state en- 
The d+Au results 



shed light on these nontrivial effects and are necessary to 
make a firm statement about thermal photon emission 
in Au+Au collisions. The CNM effects were studied in 
d+Au collisions at these energies through measurements 
of 7r°,?7 and J/ip [l2-14|; however, direct photons allow 
studying the initial-state nuclear effects - without the 
ambiguities of the hadronization process. 

In this paper, we present results of direct-photon mea- 
surements in y/s NN =200 GeV d+Au collisions at midra- 
pidity for 1 < px < 16 GeV/c. Both virtual-photon 
and real-photon measurements are performed as inde- 
pendent analyses. The virtual-photon analysis uses data 
taken in 2008 to provide results for the low pr region, 
approximately 1 < px < 6 GeV/c. The real-photon 
analysis uses data recorded in 2003 for complimentary 
results above 5 GeV/c. In addition, we report improved 
direct photon results in ^fs =200 GeV p+p collisions for 

1 < pr < 5 GeV/c using 2006 data. The new p+p results 
are combined with the previously published p+p collision 
data from 2005 to serve as a reference for the d+Au 
data. 

The two central arms of the PHENIX detector (l5| 
cover \rj\ < 0.35 in pseudorapidity and tt/2 in azimuthal 
angle for each arm. Minimum bias (MB) events were trig- 
gered by beam-beam counters located at both sides of the 
interaction point, covering 3.0 < |?7| < 3.9, which were 
also used to determine the event centrality for d+Au col- 
lisions. Events containing high px photons and electrons 
were selectively recorded by photon and single electron 
triggers in coincidence with the MB trigger. The photon 
trigger required an energy deposition in the electromag- 
netic calorimeter (EMCal) and the electron trigger re- 
quired a hit in the ring imaging Cerenkov detector with 
a correlated, above threshold, EMCal energy deposition. 
The virtual-photon analysis used 0.7 nb _1 of MB data 
and 54.9 nb _1 of single-electron-triggered data. The an- 
alyzed MB and singlc-photon-triggered data samples for 
the real-photon analysis were 0.8 and 1.6 nb _1 , respec- 
tively, where 1 nb _1 of d+Au collisions corresponds to 

2 x 197 nb _1 of nucleon-nucleon collisions. We also ana- 
lyzed 4.0 pb _1 of the p+p data from the 2006 run to mea- 
sure the direct photon cross section for 1 < pr < 5 GeV/c 
through the virtual photon analysis. 

Electron tracks above 0.2 GeV/c momentum are re- 
constructed using drift and pad chambers in each of 
the central arms, with momentum resolution a PT /px = 
1.1% © 1.16% x pr- Electrons are identified by requiring 
hits in the ring imaging Cerenkov detector and matching 
the momentum with the energy measured in the EMCal. 
Electron pairs are used to measure virtual photons using 
the method described in Ref. [H, @] . 

Any source of real direct photons also produces nearly- 
real virtual photons, i.e. low mass e + e~ pairs, allowing 
extraction of the real direct photon yield from low mass 



e + e _ pairs. In the virtual photon analysis, e + e~ pairs 
with m ee < 0.3 GeV/c 2 and pair px > 1 GeV/c are 
measured by the two central arms. Electron pairs are 
formed from combinations of all electrons and positrons 
with pt > 0.3 GeV/c in an event, and background pairs 
arising from random combinations, external conversions, 
correlated background from double Dalitz decays of 7r°, 77 
and jet induced correlations are removed by analysis tech- 
niques as discussed in Ref. Q. Electron pair mass dis- 
tributions for different pair px ranges, which comprise 
the virtual direct photon signal and the hadron decay 
component, are obtained. The inclusive photon yield 
is determined from the yield of e + e~ pairs in m ee ~ 
0.05 GeV/c 2 with the relation of = f^-^-dn 7 @j. 

The e + e~ mass distribution for m ee < 0.3 GeV/c 2 and 
Pt > 1 GeV/c is decomposed by a two-component fitting 
procedure described in Ref. Q using the known shapes 
of the direct photon and hadron decay components. The 
direct photon fraction, r 7 = direct 7/inclusive 7, is ex- 
tracted from the fitting. Multiplying the direct photon 
fraction by the inclusive photon yield leads to the direct 
photon yield. 

The systematic uncertainties on the direct-photon frac- 
tion are estimated from the difference in extracted direct- 
photon fraction when varying: (1) the particle composi- 
tions in the "cocktail" of hadron decay contributions for 
the fit, (2) the background subtraction of the measured 
mass distribution, (3) the mass region used for the fit, 
and (4) the efficiency corrections. The largest uncertainty 
is due to the particle composition of the hadronic cock- 
tail, particularly 77/71-° = 0.48 ± 0.03 atp T > 2 GeV/c, 
which is essentially identical to p+p [16|. The result- 
ing uncertainty in the direct-photon fraction due to 77/71"° 
is about 20-30%, and less than 5% are from all other 
sources. The uncertainty in the e + e~ pair acceptance 
correction introduces an additional 9% uncertainty to the 
inclusive photon yield, which is added in quadrature with 
the other uncertainties. 

Figure [T] shows the measured direct-photon fractions 
by the virtual-photon analysis in p+p, d+Au, Au+Au [l[ 
collisions from left to right. The p+p result is the com- 
bination of [l[ and the 2006 data. The curves show the 
expectations from a next-to-leading-order perturbative- 
quantum-chromodynamics (NLO pQCD) calculation [17, 
Il8j . The cutoff mass scale dependence of the calcula- 
tion is also shown for three cases: /1 = Q.hpx, 1-0 pt 
and 2.0 pr- The expectation for d+Au is calculated 
by scaling with the nuclear overlap function calculated 
from a Glauber model (l9j . which is expressed as X^a = 
^Vcoii/o-pp 6 '- Here, N co \\ is the number of binary nucleon- 
nucleon collisions and a™ el is the cross section of inelas- 
tic p+p collisions of 42 mb. The p+p data points were 
much improved statistically compared to the previously 
published data, especially above 3 GeV/c, and the p+p 
result is in good agreement with the NLO pQCD expec- 



5 



(a) p+p 



NLO pQCD 
_(i = 1.0p T 
— H = 0.5p T 
H = 2.0p T 





P T (GeV/c) 



p T (GeV/c) 



2 3 4 5 6 

P T (GeV/c) 



FIG. 1: (color online) The direct-photon fractions from the 
virtual-photon analysis as a function of pr in (a) p+p, (b) 
d+Au, and (c) Au+Au (MB) 0] collisions. The statistical and 
systematic uncertainties are shown by the bars and bands, re- 
spectively. The curves show expectations from a NLO pQCD 
calculation [l?], [3 with different cutoff mass scales: (solid) 
/i = 0.5 pr, (dash) fj, = f .0 pr, and (dash-dot) fi = 2.0 pr- 



tations. The c?+Au data are higher than the expectation 
in pt < 4 GeV/c. Their px dependence is similar to the 
NLO pQCD expectation, unlike the Au+Au data. 

While the virtual-photon analysis suffers from low 
statistics at px > 5.0 GeV/c, the real-photon analy- 
sis is robust at high px- The primary detector for the 
real-photon analysis is the EMCal, which comprises six 
sectors of lead-scintillator calorimeter and two sectors 
of lead-glass calorimeter. Contamination from charged 
hadrons is eliminated by a track-matching veto in the 
drift chamber as well as a profile cut on the EMCal 
shower. Analysis details have been described in [3I [2fjj. 
The key to the method is the precise subtraction of the 
large photonic background originating from hadronic de- 
cays, about 80% of which come from tt° — > 2j and about 
15% from r\ — > 2-f. Two techniques, 7r°-tagging and sta- 
tistical subtraction methods, are used to remove decay 
photons. 

The 7r°-tagging method identifies neutral pions by re- 
constructing pairs of photons in the lead-scintillator EM- 
Cal sectors that deposit more than 150 MeV. All pairs 
of photons at least 10 towers («0.1 radian) inside the 
edge of the EMCal which reconstruct to invariant mass 
105 < m 77 < 165 MeV are tagged as n° decays. The 
number of direct photons, jdir, is determined as 



Idv 



: lincl — (1 + -R/i/7r°)(l + ^m-iss)77r ->2 7 ! 



(1) 



where ji nc i, 7 7r o_ ! .27 are the number of inclusive and tt° 
decay photons, respectively, and Rh/-K° is the ratio of 
other hadronic contributions to 7r° decay photons. 5 m iss 
represents the probability that either of the photons from 
7T° — > 27 misses the detector. A fast Monte Carlo (MC) 
simulation, which includes the geometric acceptance and 



EMCal response, is used to estimate S m i SS . The input px 
distribution of 7r° is taken from p+p collisions 21]. 5 m i SS 
is then determined as a function of px and its uncertainty 
is evaluated as ~6-8% by varying the implemented sim- 
ulation conditions. Rh/w° is calculated using the yield 



ratios of 77 and uj to n° measured by PHENIX [2l|, |22 1 . 

The statistical subtraction method @, [23[ is applied 
to MB triggered data from both the lead-scintillator and 
lead-glass EMCal. The hadron decay contribution is es- 
timated by a hadronic cocktail simulation based on the 
observed px spectrum of 7r°; other particle spectra are 
based on the tt° using mx scaling [6|. The acceptance 
and shower merging effects are also implemented in the 
simulation. A double ratio, i? 7 , is calculated as 



R~ 



( dN-Jdpx 
ydN^o^y/dp 



data 



dNj/dpj 



dN^ 



>2~< 



jdpx 



(2) 

An excess due to direct photons gives i? 7 > 1 , and the di- 
rect photon yield is determined by jdir = (1 — R^^linci- 
Figure [5] shows the direct photon cross sections in p+p 
and d+Au collisions from both virtual- and real-photon 
analyses [4j. The NLO pQCD calculations agree with the 
p+p data well for a wide px range, and show a preference 
for the choice (x = 0.5px- Unfortunately, the NLO pQCD 
calculation with a low mass cutoff scale less than 1.0 px 
is not available for px < 2.0 GeV/c. Thus, we use an 
empirical parameterization, Eq. [3j inspired by a NLO 
pQCD formulation for p+p -J> jX fli| : 



E 



dp 3 



■Px 



(6+c-ln xt) 



(1-4)' 



(3) 



where a, b, c, and n are free parameters and xx = 
2pxj\fs- The first factor, p T < ' b+c lnXT \ is a power law 
with a logarithmic scaling correction. The convolution 
of two PDFs in colliding protons consequently introduces 
the factor, (1 — x T ) n , which naturally leads to a drop of 
the cross section to at it = 1. The virtual-photon 
(1.5 < px < 5 GeV/c) and real- photon (px > 5 GeV/c) 
results are fit simultaneously, and the point-to-point un- 
certainty of the data is considered at fitting. The px- 
correlated uncertainty of the fit is identical with that of 
the data. The quadratic sum of these fit uncertainties is 
indicated as dotted lines in Fig. [2] The fit describes the 
data very well for the entire px range. The fit param- 
eters with uncertainty (excluding the pT-correlated un- 
certainty) are a=6.6±3.3) x 10"3, 6=6.4±0.3, c=0.4±0.2, 
and n=17.6±14.9, with x 2 /NDF=22.4/16. The factor 
of the power law, b + c ■ lnix, becomes 4.6-5.5 for 
0.01 < x T < 0.1. 

The d+Au data illustrate full consistency between the 
three aforementioned independent analyses. The inde- 
pendent results are in good agreement in the overlap re- 
gion from 3.0 < px < 6.0 GeV/c. The virtual photon 
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FIG. 2: (color online) fa) The invariant cross sections of the 
direct photon in p+p @, Q arm d+Au collisions. The p+p fit 
result with the empirical parameterization described in the 
text is shown as well as NLO pQCD calculations, and the 
scaled p+p fit is compared with the d+Au data. The closed 
and open symbols show the results from the virtual photon 
and 7r°-tagging methods, respectively. The asterisk symbols 
show the result from the statistical subtraction method for 
d+Au data, overlapping with the virtual photon result in 
3 < pr < 5 GeV/c. The bars and bands represent the 
point-to-point and p-r-correlated uncertainties, respectively, 
(b) The p+p data over the fit. The uncertainties of the fit 
due to both point-to-point and pr-correlated uncertainties of 
the data are summed quadratically, and the sum is shown as 
dotted lines. The NLO pQCD calculations divided by the fit 
are also shown. 



analysis reaches down to 1 GeV/c, and the 7r°-tagging 
method extends to 16 GeV/c. The d+Au data are in 
agreement with the binary collision scaled p+p fit result 
across the entire pt coverage. A power law fit, Ap^, n , 
is performed with the d+Au data for pt > 8 GeV/c as 
done for p+p (n = 7.08±0.09 stat ±0.1 syst ) H and Au+Au 
(n = 7.18±0.14 stat ±0.06 syst for most central) @. The fit 
gives a power of n = 7.17 ± 0.76 stat ± 0.01 syst , consistent 
with p+p and Au+Au. 

Figure [3] shows the nuclear modification factor for 
d+Au, RdA, calculated as the d+Au data divided by 
the binary-scaled p+p fit. The point-to-point and pr- 
correlated uncertainties of the p+p fit are quadratically 
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FIG. 3: (color online) Nuclear modification factor for d+Au, 
RdA, as a function of pr- The closed and open symbols show 
the results from the virtual- and real-photon measurements, 
respectively. The bars and bands represent the point-to-point 
and pr-correlated uncertainties, respectively. The box on the 
right shows the uncertainty of T^a for d+Au. The curves 
indicate the theoretical calculations [24| with different combi- 
nations of the CNM effects such as the Cronin enhancement, 
isospin effect, nuclear shadowing and initial state energy loss. 



summed with those of the d+Au data points. The sums 
are shown as bars and bands, respectively. The uncer- 
tainty of TdA for d+Au is indicated by the box located at 
the right in Fig. [3] RdA is consistent with unity within 
the reported uncertainty. The theory calculations [24| 
with different combinations of standard CNM effects are 
shown with the data. The solid curve is the simplest 
one including only the Cronin enhancement and isospin 
effect. The nuclear shadowing with the EKS98 parame- 
terization 25[ of the nPDFs is additionally considered for 
the dotted and dash-dotted curves, and the initial state 
energy loss is included for the dash-dotted curve. The 



data are consistent within uncertainties with the theo- 
retical calculations, but do not have a sufficient preci- 
sion to resolve the considered initial state nuclear effects. 
The data do however rule out much larger effects be- 
yond these standard range predictions. In contrast, Fig.0] 
shows that for Raa in Au+Au collisions, there is a much 
larger enhancement of the direct photon production be- 
low 2.0 GeV/c. The magnitude of the enhancement in 
Au+Au with Raa > 7 is much higher than observed in 
d+Au, indicating that there is a significant medium effect 
on direct photon production. 

In conclusion, direct photons in 1 < pt < 16 GeV/c 
have been measured for d+Au collisions via three in- 
dependent analyses, the virtual photon, 7r°-tagging and 
statistical subtraction methods. The results from these 
analyses agree in the overlap px region. The p+p spec- 
trum has also been improved statistically by the 2006 
data. The improved p+p data are parameterized by a 
pQCD inspired fit function. The fit describes the data 
very well for the entire px region. RdA is consistent with 
unity. The data fully support the theoretical calculations 
with the standard CNM effects for a wide pr range. Raa 
shows a much larger enhancement below 2.0 GeV/c com- 
pared to the d+Au data, indicating the existence of a 
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FIG. 4: (color online) Nuclear modification factors for Au+Au 
(MB) and d+Au as a function of pr- The triangle symbols 
show results from the (closed) virtual [l| and (open) real pho- 
ton [f§| measurements, respectively. The bars, bands, and box 
represent the same uncertainties as in Fig. [3] The (+) sym- 
bols for RdA for pr < 5 GeV/c illustrates the difference in 
magnitude for -Raa between Au+Au and d+Au collisions. 

medium effect as an additional source of direct photons. 
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